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Current cosmological observations, when interpreted within the framework of a homogeneous and 
isotropic Friedmann-Lemaitre- Robertson- Walker (FLRW) model, strongly suggest that the Universe 
is entering a period of accelerating expansion. This is often taken to mean that the expansion of 
space itself is accelerating. In a general spacetime, however, this is not necessarily true. We attempt 
to clarify this point by considering a handful of local and non-local measures of acceleration in 
a variety of inhomogeneous cosmological models. Each of the chosen measures corresponds to 
a theoretical or observational procedure that has previously been used to study acceleration in 
cosmology, and all measures reduce to the same quantity in the limit of exact spatial homogeneity 
and isotropy. In statistically homogeneous and isotropic spacetimes, we find that the acceleration 
inferred from observations of the distance-redshift relation is closely related to the acceleration of 
the spatially averaged universe, but does not necessarily bear any resemblance to the average of the 
local acceleration of spacetime itself. For inhomogeneous spacetimes that do not display statistical 
homogeneity and isotropy, however, we find little correlation between acceleration inferred from 
observations and the acceleration of the averaged spacetime. This shows that observations made in 
an inhomogeneous universe can imply acceleration without the existence of dark energy. 



PACS numbers: 98.80.Jk 



I. INTRODUCTION 



The question of if, and why, the expansion of the 
Universe is accelerating is one of the foremost prob- 
lems in fundamental physics today. Observations 
of distant Type la supernovae are well-fit by homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) solutions of general relativity, but only 
if a substantial fraction of the energy density in the Uni- 
verse is in the form of a cosmological constant. For 
this constant to take the value required by observations, 
however, requires an extraordinary degree of fine-tuning. 
Many authors have attempted to remedy this situation 
by proposing modifications to gravity, the existence of 
negative-pressure components of the cosmological fluid, 
new scalar fields, and a host of other exotic mechanisms. 
Others have sought to explain it by modifying the basic 
assumptions of cosmology itself, such as the Cosmological 
Principle. 

Our aim here is to clarify what is meant by acceleration 
in an inhomogeneous universe, and to study how differ- 
ent types of acceleration can arise within the context of 
relativistic cosmology. We do not consider the existence 
of any strange new matter fields or modifications to Ein- 
stein's equations, but we do allow for the Universe to be 
strongly inhomogeneous below a certain scale. It is found 
that observations made over large distances (e.g. using 
supernova and CMB observations) are best modeled us- 
ing non-local averages of geometric quantities. The accel- 
eration that one can infer from these large-scale probes 
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of cosmology is not necessarily similar to any locally de- 
fined measures of acceleration, and does not necessarily 
imply that the expansion of space itself is accelerating at 
any point. 

Following Hirata & Seljak [1] and Clarkson & Umeh [2], 
we consider multiple different measures of acceleration: 

(a) . The local acceleration of nearby observers in a small 
volume of space, as governed by the Raychaudhuri 
equation; 

(b). The acceleration inferred by fitting an FLRW 
distance-redshift relation to Hubble diagrams con- 
structed from observations made over large dis- 
tances; 

(c) . The acceleration inferred by constructing a Hubble 
diagram from local observations (i.e. at z ~ 0), as 
described by Kristian and Sachs ^SJ; 

(d) . The average acceleration of a large volume of space, 
as calculated using Buchert's scalar averaging pro- 
cedure [4]. 

For an exactly homogeneous and isotropic FLRW model, 
all of these measures are identical. In the general case, 
however, they can be quite different, to the point that 
some of them can show strong acceleration, while others 
show deceleration. 

Measures (a) and (c) are purely local, and depend only 
on the curvature of spacetime within the neighbourhood 
of a point (in the limit of geometric optics). Measures 
(b) and (d), however, are both inherently non-local in 
nature. Measure (a) describes the actual, dynamical ac- 
celeration of a volume of space, while measure (b) corre- 
sponds to what is measured in, for example, supernova 
surveys. Measure (c) is a local measure of acceleration. 



based on inferences made using the luminosity distance- 
redshift relation within a very small region of spacetime. 
Measure (d) corresponds to the theoretical procedure of 
fitting a smooth effective model to the real, inhomoge- 
neous Universe (the question of how best to do this is still 
very much an open one, but the simplest and best-known 
procedure is the scalar averaging formalism developed by 
Buchert OH]). We will discuss all of these measures in 
detail in what follows. 

Constructing relativistic models of the real, inhonio- 
geneous Universe is difficult, but a number of attempts 
have met partial success in addressing the problem. Re- 
alistic density fields can be described using perturbations 
to exact FLRW solutions. However, at least to quadratic 
order in perturbations, the effect of inhomogeneities on 
the overall behaviour of the spacetime often turns out to 
be small, with the average evolution and optical prop- 
erties of the model remaining close to the unperturbed 
background values pilTHTT]. There have been suggestions 
that effects at higher orders in perturbations could be im- 
portant f2| [12] , and recent studies of non-linear collapse 
which make use of a gradient expansion also find poten- 
tially significant effects [IHl US] . These claims require fur- 
ther analysis if they are to be confirmed, however. The 
main drawback of the perturbative approach is that the 
solutions obtained are not exact, and so it is unclear if 
we are neglecting or incorrectly estimating aspects of the 
fully-relativistic behaviour of the spacetime. 

An alternative approach involves studying exact, in- 
homogeneous solutions to Einstein's equations. As the 
geometry of spacetime is known from the outset, the non- 
linear evolution of space, and light rays within it, can be 
calculated in a fully rigorous manner, without recourse 
to perturbative analyses. This approach tends to require 
a high degree of symmetry, however, placing strong re- 
strictions on what it is possible to model [2^. For ex- 
ample, backreaction effects in Swiss Cheese models (con- 
structed by joining together spherically-symmetric inho- 
mogeneous regions with a FLRW solution) are known to 
be small [21]. This is confirmed by ray tracing studies 
through a variety of similar exact solutions pT| - B5] . Ad- 
ditionally, inhomogeneous models with general fluid con- 
tent are difficult to work with [251 - 151] , and so most mod- 
eling attempts make do with a dust-only stress-energy 
tensor. This is itself problematic, as high-density regions 
tend to rapidly collapse in inhomogeneous pressureless 
fluids, forming singularities and resulting in unrealistic, 
pathological behaviour. 

Finally, more heuristic analyses have been attempted 
[3TI [35H35] . These often employ disjoint regions of differ- 
ent exact solutions. The dynamical evolution within each 
region is therefore well deflned, but the ensemble as a 
whole does not satisfy Einstein's equations at the bound- 
aries between regions. This approach relaxes the restric- 
tive symmetry requirements that are necessary when us- 
ing exact solutions, but introduces its own ambiguities. 
Indeed, the chosen boundary conditions can sometimes 
dominate the behaviour of the model 1331. and it is not 



clear whether the resulting spacetime approximates any 
actual solution of Einstein's equations or not. 

In this paper we consider the latter two of these ap- 
proaches, as the situation of small perturbations around 
an exact FLRW geometry is well-studied already O IT6l 
IHSHHj- The format of the paper is as follows. In Section 
im we discuss in detail the different measures of accel- 
eration summarised above, and how they can be calcu- 
lated in a general spacetime. The conditions for these 
measures to show acceleration, and their relation to ob- 
servable quantities (if any), is discussed. Section III sets 



out three different inhomogeneous models: the spherical 
collapse model, constructed from disjoint FLRW regions; 
the Kasner-EdS model, an exact solution with alternat- 
ing expanding vacuum and collapsing dust regions along 
a line of sight [431 - 147] : and the Lemaitre-Tolman-Bondi 
model, an exact, spherically-symmetric dust solution [481 - 
[50] . The volume of space in all of these models is locally 
decelerating everywhere, and yet we find that observa- 
tions made within them can still exhibit acceleration. 



In Sections |IV[ \V\ and [VI[ we present our results 
for each of the models. We find that the acceleration 
inferred by observers from the Hubble diagram, con- 
structed over large distances, is most closely correlated 
with the dynamical behaviour of the averaged spacetime, 
and not with its local acceleration. The local accelera- 
tion, however, is well modeled by Hubble diagrams that 
are constructed locally, using the Kristian-Sachs formal- 
ism. This result suggests that we should use non-local 
averages of the geometry to interpret observations made 
over large distances, and local geometry to interpret ob- 
servations made locally. Interpreting observations within 
the wrong framework can lead to incorrect inferences 
about what is causing the acceleration. This appears to 
us to provide some insight into the question of whether 
the apparent acceleration we observe is necessarily caused 
by the presence of a non-zero cosmological constant, or 
whether it could be caused by some backreaction effect 
within the inhomogeneous spacetime. Finally, in Sec- 
tion [VlT] we discuss our results in the context of previous 
claims involving backreaction and the nature of the ap- 
parent acceleration of the Universe. 



II. MEASURES OF ACCELERATION 



In this section we discuss the four measures of accel- 
eration that were listed in Section [l] Each measure has 
associated with it a deceleration parameter. For each def- 
inition, we set out its theoretical basis and physical inter- 
pretation, list the conditions necessary for acceleration to 
occur, and describe the relation of the deceleration pa- 
rameter associated with this measure to observable quan- 
tities. The measures are summarised in Table |T] at the 
end of the section. 



(a) Local volume acceleration 

A congruence of time-like geodesies, u°-, describing the 
world-lines of a set of observers comoving within a cos- 
mological fluid can be decomposed such that j51j 



Ua;b 



l^'dhab 



- (^ab +^ab, 



(1) 



where subscript ; denotes a covariant derivative, and 
hab = gab + UaUh is the projection tensor. The kinematic 
quantities in this equation are the expansion scalar, 0, 
the shear, aab, and the vorticity, tOab, which correspond to 
the trace, symmetric trace-free, and antisymmetric parts 
of Ma;f), respectively. This decomposition is a fully covari- 
ant procedure, valid for any spacetimc. 

For irrotational flows {ojab = 0) , the fluid flow becomes 
hypersurface orthogonal, with the projection tensor be- 
coming the induced metric of the orthogonal 3-spaces. 
The field equations then give (with 8ttG = 1) 



e 



-^e2-i(p + 3p) + A-2a2, (2) 



where the over-dot denotes a derivative with respect to 
proper time along u°, and where a^ = ^CTabcr"^, and p, p 
and A are the energy density, pressure and cosmological 
constant, respectively. The Gauss embedding equation 
for these hypersurfaces is 



1 



e^=p + A 



1 



(3)-i 



'y'-^'"n, 



(3) 



where *^^^7?. is the Ricci scalar constructed from hab- 

By analogy with the deceleration parameter in an 
FLRW spacetime, one can then define a local volume 
deceleration parameter. 



qe 
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-{p + 3p)-A + 2a' 



(4) 



corresponding to the monopole of the deceleration of the 
expansion rate of a set of neighbouring particles following 
m". This measure is said to be accelerating when qq < 0, 
which from Eq. Q can be seen to occur if and only if 
p < -|p or A > 0. 

Now, qq is a local measure of acceleration, defined only 
in the neighbourhood of a point in spacetime. Deter- 
mining ge from observations therefore requires (and 
its first derivative) to be determined using observations 
within the neighbourhood of a single point in spacetime 
only. These observations need to be of the rate of change 
of proper distance between particles in their own rest- 
frame, which is not necessarily the same as the angular 
diameter distance or luminosity distance (for these see 
measure (c), below). The Universe is extremely inhomo- 
geneous on small scales, and so this measure of accelera- 
tion is likely to display considerable spatial variation. A 
local measurement of qq need not, then, be representa- 
tive of the mean local volume acceleration. 



(b) Observed acceleration 
(as inferred from the Hubble diagram) 

To date, the strongest observational evidence for an ac- 
celerating universe comes from the distance-redshift re- 
lation, which is measured using "standardisable candles" 
(such as Type la supernovae), as well as the CMB. By 
fitting this data to the relations derived from the FLRW 
solutions of Einstein's equations, one finds that models 
with 51a > are strongly favoured, while those with 
riA = are ruled out to a high degree of confidence. Un- 
der the assumption that spacetime is well-described by 
FLRW, this constitutes strong evidence for accelerating 
expansion. 

A key step in this procedure is the fitting of the FLRW 
distance-redshift relation to data. Performing a series 
expansion in the angular diameter distance about z = 
gives 



dA{z) 



cz 



--{■i + qo)z + 0{z') 



(5) 



where Hq is the inferred Hubble rate at z = 0, and go 
is the inferred deceleration parameter. The value of qo 
can be found in a purely observational manner by taking 
derivatives of the fitted distance-redshift relation W^ , 



Qo 



rl" 

d' 



(6) 



where primes denote derivatives with respect to redshift, 
and where subscript denotes that a quantity is evalu- 
ated at the observer. Note that this measure of deceler- 
ation is inferred from observations over large distances, 
and so assigns information to a point in spacetime based 
on information obtained from the entire extended region 
over which observations have been made. This is a highly 
non-local process. 

Eq. (pi) is derived from the distance-redshift relation 
in FLRW geometry, but one should note that the exis- 
tence of such a geometry is not required in order to infer 
go from observations. What we are doing here should in- 
stead be considered simply as a fitting procedure. Eq. ^ 
is then an observable in any spacetime (after averaging 
over the celestial sphere), and so anyone can measure go 
if they are willing to interpret their observations within 
the framework of an FLRW model. We call the accel- 
eration inferred by this fitting of FLRW relations to the 
monopole of the distance-redshift relation the observed 
acceleration. 

The distance-redshift relation can be calculated in a 
general spacetime by solving the Sachs optical equations 
for bundles of null geodesies. With vanishing vorticity, 
these are 



de 

dX 



1 



aub 



da 
dX 



2ea 



Jriabf^ fc 



Cabcd{t*rk''{t*rk'', 



(7) 
(8) 



where A is an afRne parameter along the bundle, Rab 
and Cabcd are the Ricci and Weyl tensors of the space- 
time, 9 and a are the expansion and (complex) shear 
scalars of the null geodesies, k"' is a tangent vector to 
the null curves, and i" are (complex) vectors spanning 
a two-dimensional screen space orthogonal to k°'. The 
expansion scalar is related to the angular diameter dis- 
tance measured along the bundle hy 9 = d{\n{dA))/dX. 
Substituting this into Eq. ^ yields 



dA2 



^dA 



lap 



1 



JrCabf^ 1^ 



(9) 



The luminosity distance is related to the angular diame- 
ter distance by the reciprocity theorem [53 , which gives 
dL = {l + zfdA. 

The afhne distance-redshift relation for a general 
spacetime is given by j54j 



dz 
dX 



-(l + z)2i/||(z), 



(10) 



where iJ|| = |8-f o'abe^e'' is the expansion rate along the 
line of sight, e°. Equations ([8|-([Iol) can be solved to give 
the angular diameter distance-redshift relation, (i^(z), in 
any given direction, at any given point in spacetime. This 
procedure can be repeated for every direction on the sky, 
and an FLRW model can be fitted to the monopole of 
the resulting angular distribution. The best-fit model 
can then be used to find go • We denote the result of this 
procedure in a general spacetime as 9obs- 

We consider this definition of acceleration to corre- 
spond most closely to the one used by observers. It is 
a non-local measure, since it depends on solutions to the 
Sachs equations, which describe bundles of null geodesic 
curves that extend through the spacetime. In effect, 
the measurement of gobs depends on finding the entire 
past null cone of an observer out to some z, and fitting 
some distance-redshift relation to it. The conditions for a 
spacetime to have gobs < are therefore complicated, in 
general. As we will see below, it is possible to find space- 
times that are quite different from simple A-dominated 
FLRW models that nevertheless have gobs < 0. 



(c) Acceleration from local observations 
(using the Kristian-Sachs formalism) 

The measure of acceleration we just described has the 
disadvantage of requiring solutions to the Sachs equations 
to be found, as a function of redshift and angle on the 
observer's sky. This can be a difficult task in general, as 
it requires detailed knowledge of the geometry of space- 
time. Instead, one can use the Kristian-Sachs formalism 
[3] to obtain a fully general and covariant series expansion 
of the distance-redshift relation about an observer with- 
out needing to consider solutions to the geodesic equa- 
tions at all. Furthermore, the expansion can be decom- 
posed directly into covariant spherical harmonics about 



the observer, allowing the monopole term to be calcu- 
lated straight away [2]- 

This procedure has been spelled out in detail by Clark- 
son and Umeh [2j . The generalised form of Eq. (Is]) is 



dA 



[K-KbWaUb]o 



K'^K^K'VaVbUc 



(11) 
where the past-pointing null direction can be written in 

terms of the tangent vector to a comoving observer's 

world-line, m°, and a direction on their sky, e°, as 



K" = 



k'' 



[ubk% 



= -u" + e" 



(12) 



The subscript again denotes evaluation at the ob- 
server's location. The terms in Eq. ( 11 ) can be expanded 



using a covariant decomposition in spherical harmonics. 
In order to facilitate this expansion, it is useful to invert 



Eq. (11) to give 



1, 



[K''K''VaUb]odA + 2 [i^^i^'if "VaVfcwJod; 



i\) 



0{d', 

. (^^) 
All of the terms that we wish to expand are now in the 

numerator. Comparing the monopoles of the coefficients 

in Eq. ( 13 ) with the corresponding FLRW relation (with 

SttG = 1) then gives 



gKS 



3 

62 



(p + 3p) - A + Qa' 



(14) 



where we have used Hq = l9, which corresponds to 
the monopole of this term, rather than its full spherical 
harmonic expansion. This corresponds most closely with 
the way that iJo is typically used in observational studies; 
the monopole of Hq tends to be determined separately 
from other quantities. From Eq. ( 14 ), it can be seen that 
gKS > unless p < — |p or A > 0. That is, acceleration 
of this measure can only occur if there is a cosmological 
constant, or if an exotic fluid with negative pressure is 
present. 

There is clearly some similarity between the decelera- 
tion parameter of the local volume, ge, and that which 
is obtained from the local distance-redshift relation, gxs- 
Both measures of acceleration are local (depending only 
on quantities defined within the neighbourhood of the 
observer), and both are given by expressions that dif- 
fer only by a term involving the shear scalar. This does 
not, however, mean that these two measures of accelera- 
tion are the same. They correspond to different physical 
quantities. 

The relation between the Kristian-Sachs and observa- 
tional acceleration measures is less clear. This has to do 
with the non-locality of the latter; if the Hubble diagram 
could be measured precisely at z = (i.e. in the limit of 
geodesies of zero affine length), we would find gobs = gxs- 
But, because real observations necessarily cover a range 
of redshifts, the process of fitting a curve to the data and 



extrapolating that back to z = means that in general 
we will have gobs 7^ ^ks- This phenomenon has been in- 
vestigated using real data in [55] • As we shall see in what 
follows, even the signs of gobs and gxs can be different. 
That is, a locally-decelerating spacetime can still have a 
Hubble diagram that implies acceleration. 



(d) Acceleration of the average 
(using Buchert's formalism) 

One method of constructing a homogeneous and 
isotropic 'effective' model within which observations can 
be interpreted involves taking averages of geometrical 
quantities over space- like hypersurfaces. The hope is 
then that the behaviour of the averaged model will cap- 
ture some aspects of the real spacetime, both in terms of 
its dynamics, and the observational quantities that are 
calculated within it. Many such averaging procedures 
exist [56 , but here we will concentrate on the Buchert 
scalar averaging formalism 4_. This is the most widely 
used formalism in the literature. 

Buchert's method proceeds as follows. First of all, 
the spacetime is filled with a congruence of (irrotational) 
curves. It is then foliated with a set of space-like hyper- 
surfaces orthogonal to these curves. The proper 3- volume 
of a domain, P, on a given hypersurfaces is 



Vn 



■D 



(15) 



where hij is the induced metric on the hypersurface, and 
h = det hij. In general, the induced metric will be a 
function of time, and so the volume is time-dependent 
as well. The proper volume-weighted average of a scalar 
quantity, S*, over a spatial domain, P, can then be written 
as 



{S) = V, 



V 



S{x,t)V^ d^x. 



(16) 



V 



Spatial averaging and time evolution do not, in general, 
commute. They instead obey the commutation relation 



dt{S)-{dtS) = {QS)-{Q){S). 



(17) 



An 'effective' homogeneous model can be constructed by 
averaging over domain sizes greater than the statisti- 
cal homogeneity scale of the underlying inhomogeneous 
spacetime. An effective scale factor for the resulting 
model can then be defined as 



av{t) 



( VvJt) 
\Vv{to) 



(18) 



where Iq is some fiducial time. For geodesic curves, and 
pressure-free matter, we can use this formalism to con- 
struct analogues to the Friedmann and Raychaudhuri 



equations, 
3^ 



8^G(p) +K--[Qt> 
-87:G{p)+A + Qt,, 



/(3) 



7^) 



(19) 



where -//d = dx>/ax>. Over-dots denote partial differen- 
tiation with respect to t, the proper time along curves 
orthogonal to the 3-space. The kinematical backreaction 
scalar is defined to be 



{{Q')-{ef)-2{<j'). 



(20) 



A deceleration parameter for the averaged hypersurfaces 
can then be defined, by analogy with FLRW cosmological 
models, as 



Qv 



1 ap 
H^ aj,' 



(21) 



An effective distance-redshift relation can also be found 
for the averaged model by assuming that light rays fol- 
low null geodesies of the averaged spacetime, and that 
geodesic observers are comoving in the average geome- 
try. 

In general, the deceleration parameter qx> is non-local, 
and not directly observable, as it depends on averages 
over extended space-like hypersurfaces. From Eq. (19 1, 



one can see that if the averaged spatial Ricci curvature, 
{'^^'TZ), and backreaction, Q-p, behave in a certain way, 
then it is possible to have q-p < without having A > or 
(p) < 0. In particular, spacetimes consisting of collaps- 
ing regions in an expanding background can exhibit this 
behaviour, which has led some to claim that the apparent 
cosmic acceleration inferred from supernova observations 
could instead be explained as a consequence of the vari- 
ance of the inhomogeneous expansion rate that enters 
into the definition of Qx> [3 E |36l [Ml EH [58] . We will 
examine this claim in detail later on, but for now will just 
note that it is possible for the Buchert averaged model 
to accelerate, even if ge > everywhere. 



Measure 


Local 


Support 


Observable 


Local volume, qe 
Observational, gobs 
Kristian-Sachs, qks 
Buchert average, q-o 


Yes 
No 
Yes 
No 


Spacetime point 
Null geodesic 
Spacetime point 
Spatial domain 


In principle 

Yes 

In principle 

No 



Table I: Summary of the different measures of acceleration 
defined in Section Ull 

Quantities of particular interest are the deceleration 
parameters that observers in a given region of spacetime 
should expect to infer, rather than the actual values of 
these parameters at single points (which may not be rep- 
resentative). To this end, we also use Eq. (16) to aver- 
age our various measures of acceleration. 

(9KS) 



For example, 
will be taken to correspond to the value of the 



Kristian-Sachs deceleration parameter that a typical ob- 
server would expect to measure in a given region. Unless 
specified otherwise, the averaging domain is taken to be 
larger than the homogeneity scale of the model^ . 



III. INHOMOGENEOUS COSMOLOGICAL 
MODELS 

We now consider a number of different relativistic mod- 
els, and how the measures of acceleration defined in Sec- 
tion [n] can be calculated in each of them. We will pri- 
marily be interested in cases with vanishing cosmological 
constant, as a non-zero A generically results in acceler- 
ating expansion in all of the measures discussed above. 
The models we consider are chosen as illustrative exam- 
ples, as they are capable of showing acceleration in some 
measures, while displaying deceleration in others. 



A. Spherical collapse model 

In this section we consider the spherical collapse model, 
which consists of an ensemble of disjoint FLRW regions 
[551 - 15^ . It is not an exact solution to Einstein's equa- 
tions (distinct FLRW models cannot satisfy the appro- 
priate junction conditions on their shared boundary), but 
provides a useful toy model of collapsing structures in an 
otherwise expanding universe. By virtue of each indi- 
vidual region being homogeneous, spatial averages take 
a particularly simple form in spherical collapse models. 
A scalar quantity S may be averaged simply by taking 
the sum of its values in the different regions weighted by 
the proper 3-volunie of each region, V, at a given time 
0133, 



(S)^ 






(22) 



We wish to model the collapse of structures in an expand- 
ing background, as this configuration is known to lead to 
acceleration of the Buchert averaged spacetime, qv < 
[3 [31] . To this end, we select alternating regions of col- 
lapsing, spatially closed and dust-dominated regions, and 
expanding, spatially open vacuum regions, as illustrated 
in Fig. [ij We model the regions using FLRW geometry 
(with A = 0), and refer to them as Regions I and II, re- 
spectively. The collapsing regions are chosen so that they 
have a comoving size of order 10% of the vacuum regions 



It should be noted that this averaging scheme is weighted by 
proper volume, and so the density of observers is implicitly as- 
sumed to be weighted in a similar manner. Of course, this need 
not be the case, and other ways of distributing observers through- 
out space could be considered. This would affect what a 'typical' 
observer should expect to see. We do not consider this question 
further here. 



at the present time. As we look further back in time, 
the collapsing regions increasingly dominate the proper 
spatial volume. 

In what follows we will use the following form for the 
FLRW metric, 



ds' = dt' 



a2 {dX^ + dY^ + dZ^) 



(23) 



and take the individual domains to have depth Xxi- The 
proper volume of each domain is 



V, 



afdXdYdZ 



■D, 



[l + (fc,/4)(A2+y2+^2)]3' 



and the comoving depth of a domain is 

r^'^' dX 

/o 



Xz 



l + (fc,/4)(A2+r2 + ^2)- 



The effective scale factor, defined in Eq. (181, is 



av{t) 



E.^.(io) 



(24) 



(25) 



(26) 



and the Buchert average deceleration parameter qx> is 
calculated according to Eq. (21 1. 



The distance-redshift relation is required to find gobs- 
It can be found by solving the Sachs optical equations 
for a ray passing through a number of different FLRW 
regions. All that is required is to ensure continuity of 
the affine parameter. A, the redshift, z, the angular di- 
ameter distance, d^, and the expansion scalar, 0, at the 
boundary, as a ray leaves one region and enters another. 
These conditions impose no restrictions on the individual 
FLRW regions. Their size, matter content, and expan- 
sion rate can be chosen separately. We choose the cosmic 
time in individual regions such that it is continuous on 
the boundary between them. 

The Sachs equations in each region are given by 



dA2 



= -47rGp, {l + z)'^dA, 



with 4-vector tangent to the null curves given by 



k" 



l + z,±- 



l + ^fc,A2,0,0 



(27) 



(28) 



The redshift as a function of affine distance is then 



dz 
dX 



-H,{l + z)' 



(29) 



We have not used 1 + z — a^^, as the total redshift along 
a ray, z{X), is not the same as the redshift that one would 
obtain by ray-tracing through a single FLRW region. 

We propagate the light rays through each region in 
turn, until the edge of the region is reached at some co- 
ordinate distance Xi (which we call the 'length' of the 




Figure 1: Schematic representation of a line-of-sight through a universe with alternating expanding vacuum regions (clear), 
and collapsing dust regions (grey). 



region). This length is a function of time (evaluated at 
the instant the ray enters a region), chosen so that the 
ray travels a proper distance in each region proportional 
to the fraction of the proper volume taken up by that 
region on a surface of constant t, i.e. 



aiXi 
aiiXn 



Vh' 



(30) 



This choice takes into account not only the different ex- 
pansion rates of the two different regions along the line 
of sight, but also the fact that as we go forwards in time 
the number density of dust regions should be expected 
to decrease. For simplicity, we fix the (comoving) length 
of the collapsing region, x/, and then solve Eq. (301 for 
Xii on entering each new region. 

The redshift through the model is given by integrat- 
ing Eq. (29). There is a blueshift as the photons pass 



through collapsing regions, and a redshift as they pass 
through expanding regions. The distance increases de- 
spite the blueshift, and the distance-redshift relation be- 
comes multi- valued. This leads to it taking on a jagged 
appearance, although the distance remains smooth as a 
function of affine parameter. The observational deceler- 
ation parameter can be calculated by fitting an FLRW 
model with dust, curvature, and A to the jagged distance- 
redshift curve, using a simple least-squares procedure. 
The FLRW deceleration parameter, defined by Eq. (|6]), 
can then be calculated to give gobs- 

The local volume deceleration parameter in each 
FLRW region is given by Eq. (B. In our chosen spherical 
collapse model, Oa = in both regions, and Oj„ = in 
the vacuum region, so the average of the local decelera- 
tion parameter reduces to 



{IB) 



1 ^mjVl 



2Vi + Vii 



(31) 



which is positive definite (i.e. decelerating everywhere). 
The equality ge = ^ks holds, since each region is FLRW, 
and therefore has vanishing shear. 



B. Kasner-EdS model 

In the previous section we considered the spherical col- 
lapse model, which consists of disjoint regions of differ- 
ent FLRW spacetimes. This model is simple to work 
with, but only forms an approximate solution of Ein- 
stein's equations. In this section we consider an exact 



solution that is inhomogeneous along the line of sight, 
with alternating regions of collapsing dust and expand- 
ing vacuum, as illustrated in Figure [T] We take the two 
space-like directions orthogonal to the line of sight to 
span a plane symmetric subspace, and enforce statistical 
homogeneity along the line of sight only. 

In order to prevent the rapid formation of singulari- 
ties that often occurs when dust is allowed to collapse in 
general relativistic models, we will take the dust dom- 
inated regions to be locally spatially homogeneous and 
isotropic. These symmetries prevent the sudden forma- 
tion of singularities at different points in space, as every 
point is taken to be identical to every other point by fiat. 
The geometry of these regions is therefore given by the 
FLRW line-element (23), where a{t) is the scale factor 
that obeys the Friedmann equation 



SttG 



(32) 



where p oc a~^ is the energy density of the dust. In 
the vacuum regions we will take the geometry to exhibit 
translational symmetry along the line of sight, but will 
not restrict the geometry to be FLRW, as it is known 
that no solutions that satisfy the junction conditions will 
exist in this case. The geometry of this region is then 
given by 



ds-" 



-dp + bl{i)dX^ + bl(i) Uy^ + dzA , (33) 



where 6i [t) and 62 [t) are the scale factors in the directions 
tangent and normal to the line of sight, respectively. 

It can be shown that the junction conditions between 
these two regions are satisfied if we identify our hatted 
and un-hatted coordinates at the boundary, and if fc = 0, 
61 = a~^/^ and &2 = a [lH US] ■ The dust dominated re- 
gions are then spatially flat FLRW, the vacuum regions 
are Kasner, and the entire geometry is an exact solution 
of Einstein's equations [IS]- These solutions are, in fact, 
a special case of the general dust solution admitting a 
three dimensional group of space-like Killing vectors on 
two dimensional planar subspaces |231l63j . but are chosen 
such that we can have collapsing regions that do not ex- 
hibit the shell crossing singularities that tend to rapidly 
form in the general case. We can arrange for either the 
dust regions to be collapsing and the vacuum regions to 
be expanding along the line of sight, or the dust regions 
to be expanding and the vacuum regions to be collaps- 
ing along the line of sight. Here we will concentrate on 
the former, which appears to us to be more in-keeping 



with the usual picture of what is expected to happen in 
the late Universe (voids expanding, and dense regions 
collapsing). 

In order to calculate observational quantities within 
this solution, along our chosen line of sight, we will need 
to specify a set of observers. In the dust dominated re- 
gions these can be conveniently chosen to be comoving 
with the fluid, such that they follow a set of geodesic 
curves with tangent vector u° = (1, 0,0,0). In the vac- 
uum regions we will also take our observers to follow 
curves with tangent vector m° = (1,0,0,0), in the coor- 
dinates used in Eq. ( 33 ) . One should note that although 



there is no fluid in this case, so that this choice is not 
unique, it is a choice that picks out a set of curves that 
are parallel to the world-lines of observers who stay at 
the boundary between regions. They are also geodesic. 
These choices therefore correspond to a congruence of 
complete geodesic curves that fill the entire spacetime. 

Let us again refer to the dust dominated regions as Re- 
gion I and the vacuum regions as Region II. Null geodesic 
curves in these two regions, in the direction of the inho- 
mogeneity, are given by 



Cl 



fc? = ^ (a, 0,0,-1) 
fc?/ - C2 (^^, 0, 0, -a) 



(34) 
(35) 



where ci and C2 are constants, and where we have again 
taken the affine parameter that defines these tangents to 
decrease into the past. The energy of a photon following 
fc", as measured by the observers following the curves m", 
is given by 



Ej 



E 



II 



Cl 

a 
C2a 



1/2 



(36) 

(37) 



where E — —Uak""- We can also see that the Sachs optical 
equations in the two regions reduce to 



(PdAj 



= -^dA 

3a5 ' 



// 



dA2 



(38) 
(39) 



where ci is the constant from Eq. (34), which will be 



different within each individual dust dominated region. 
As always, the redshift is given by taking the ratio of 
photon energy at the time of emission and observation, 
calculated using Eqs. (36) and (37). 



The trajectories of photons can be straightforwardly 
integrated between the different regions, using Eqs. (34) 



and (35 1, and by taking the value of X{i) on leaving one 



region as its initial value on entering the next. Likewise, 
the value of E can be calculated along the null trajecto- 
ries by setting its value on entering one region as being 
equal to its value on leaving the last. This gives the value 
of the constants Ci and C2 in each of the dust and vac- 



to be integrated along the null trajectory. Integration 
of this equation again requires setting dA and 9 to be 
equal at the boundaries between regions. Following the 
prescription above, the model is uniquely specified once 
we specify three pieces of information: (i) The size of the 
vacuum regions, (ii) the size of the dust regions, and (iii) 
the time until the dust regions reach the 'big crunch'. 

Now let us consider the Buchert average of this geom- 
etry. The usual procedure is to average the expansion 
scalar over a space-like hypersurface, and use the aver- 
aged value to calculate observables. Here we have created 
a model that is inhomogeneous in one direction only. Av- 
eraging in all three spatial directions should not therefore 
be expected to reproduce anything like the observations 
we can calculate by looking in the direction of the in- 
homogeneity. Instead, we consider that the appropriate 
thing to do (and in analogy to the case of inhomogeneity 
in all 3 spatial directions), is to average the scale factor 
in the direction of inhomogeneity only. We are then left 
with an averaged geometry with line-element 

ds^ = ~dt^ + {hf dX^ + a^{t) (dY^ + dZ^) , (40) 

where the averaged scale factor, (6), is given by 



{b) = 



J y/gxxdX 

Jdx ■ 



(41) 



We can now calculate observables in this averaged geom- 
etry, and compare them to the observations made along 
the line of sight in the actual geometry of the spacetime. 
Let us now consider the spatial average of the local 
volume deceleration parameter. Again, we do not want 
to average over all spatial directions, as we are only con- 
sidering inhomogeneity and observations along one pre- 
ferred direction. The average of qq along the line of sight 
is therefore given by 



(qe) 



! qey/gxxdX 
J y/gxxdX 



(42) 



uum regions respectively, and allows Eqs. (38) and (39) 



where qq is 1/2 in the collapsing dust regions, and —4 
along the X-direction in the vacuum regions. We find a 
similar expression for gxs j which takes the same value as 
qq in both the dust and vacuum regions. 



C. Lemaitre-Tolman-Bondi model 

In the previous two sections, we considered models con- 
sisting of alternating expanding vacuum and collapsing 
dust regions. We will now consider a model with no dis- 
continuities in the density distribution, in the form of 
the spherically-symmetric, dust-only Lemaitre-Tolman- 
Bondi (LTB) solutions [48 50 . These have been the fo- 
cus of much recent interest due to their ability, in the 
guise of 'giant void' models, to reproduce the observed 
supernova Hubble diagram without a cosmological con- 
stant [B1H7D] . They are also capable of having an accel- 
erating spatial average under certain conditions [7TH7i] . 
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The LTB metric is given by 



ds' = dt-" 



aUt,r) 



(l-fc(r)r2) 



-dr^ 



al{t,rydn'^, (43) 



where 02 = {air)' is the radial scale factor. The trans- 
verse scale factor ai must satisfy the analogue of the 
Friedmann equation, 



87rGm(r) fc(r) 



A 
3"' 



(44) 



Primes and over-dots denote partial derivatives with re- 
spect to r and t, respectively. The functions k{r) and 
m(r) are arbitrary functions of the radial coordinate, and 
may be interpreted as the spatial curvature and a mass 
density at a given radius. The sign of k{r) classifies the 
differential equation, and analytic parametric solutions 
exist for each sign. Integrating Eq. ( [44| ) with respect to 
time introduces a third arbitrary radial function, ^^(r), 
which describes the local time since the big bang singu- 
larity along the world-lines of the dust. 

The metric in Eq. ( 43 ) is invariant under the trans- 
formation r -^ fir), which can therefore be used to set 
one of the arbitrary functions to a simple form, without 
losing any generality. The LTB solutions are isotropic 
about r = only, and in general have different expansion 
rates in the radial and transverse directions {Hi = hi/ai 
and H2 — 02/02, respectively). The density, expansion, 
and shear scalars for this spacetime are 



P 

e 

2 



{mr 



3V 



3a2a^r2 
2Hi + H2 
1 



;{Hi 



H,?. 



(45) 
(46) 
(47) 



Observers away from the centre of symmetry have 
anisotropic distance-redshift relations, in general. Rather 
than solving the Sachs equations in full for every direc- 
tion on the sky, we will operate within the dipole approx- 
imation^ , which assumes that the dipole term dominates 
the anisotropy of the distance-redshift relation [75ti78j . 
The dipole is aligned with the radial direction due to 
the symmetry of the model, and so the monopole can be 
estimated by taking the mean of the angular diameter 
distance in the radial directions facing into and out from 
the centre of symmetry, 



1 



dA{z)\e=o w -[dA{+r,z) + dA{-f,z)]. 



(48) 



We expect this to be a reasonable approximation for the 
models considered here. The observational deceleration 
parameter, gobs, is defined using the monopole of dA only. 



^ See the end of Section IVll and reference 
the validity of this approximation. 



TEi, for a discussion of 



For light propagation purely in the radial direction, the 
tangent vector to the null geodesies is 



k-=l{l + z),±^^^^^{l + z),0,o], (49) 



02 



and the redshift is given by integrating 



dz 
dX 



-H2{l + zf 



(50) 



The angular diameter distance in either direction can be 
found using the Sachs equation ([9]) in the radial direction. 
In general, spatial averages in an LTB spacetime will 
be both position- and domain-dependent. We consider 
spatial averaging only for spherical domains centred at 
r = 0, on hypersurfaces of constant t. We define the 
effective scale factor to be a-p = {V-d/V-d.q)^ , where 



Vn 



V 



An 



2 2 
a2afr^ 



•^1 — k{r)r'^ 



zdr 



(51) 



and rxi is the radius of the spherical domain. The average 
of a scalar quantity is then given by 



{S) = 



Att 



Vn 



V Jo 



S{r,t) 



2 2 
a2a\r^ 



^\ — k{r)r'- 



:dr. 



(52) 



We consider only the class of LTB models with A = 0. 
This means that the spacetime is locally decelerating ev- 
erywhere, ge ^ 0. The Kristian-Sachs measure is also 
necessarily decelerating (gxs ^ 0). 

Nevertheless, the freedom in the radial profiles makes 
it possible to construct models in which the Hubble di- 
agram exactly matches that of an accelerating FLRW 
model for z > 0, as seen by an observer at the centre of 
symmetry (although see [75]). As a result, LTB models 
have been studied extensively as a conventional relativis- 
tic explanation of the apparent acceleration inferred from 
supernova observations which does not require the exis- 
tence of an exotic dark energy component, or modifica- 
tions to the theory of gravity (e.g. [55 1 |64H66 l [77 l [80H85 ] ). 
Despite some success, however, they ultimately seem un- 
able to account for certain combined sets of cosmological 
observables^ [77] . 

Acceleration of the Buchert average, q-u < 0, has been 
demonstrated for a number of different LTB models (e.g. 
[7TH7?1 |5S1|57]). The cases studied generally take spher- 
ical averaging domains, centred about the origin. They 
typically find acceleration only for finite ranges of do- 
main size, which do not correspond to the homogeneity 
scale (if one exists). Additionally, models in which ob- 
servers at the centre of symmetry infer acceleration from 
their Hubble diagram seem not to correspond to those 
with an accelerating Buchert average [71] , and vice versa 



^ For an alternative perspective, see 1301 . 
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|88j . Since the spacetime is always decelerating locally, 
the existence of acceleration in the Buchert average must 
be caused purely by the backreaction term, Eq. (pO). 



In what follows, we will specialise to LTB models with 
a simple Gaussian spatial curvature profile 



k{r) = Ak exp 



(53) 



with ts ^constant, and with a choice of radial coordinate 
such that 771 (r) ^constant. When ^fc < 0, there is a cen- 
tral void region (often surrounded by an over-dense shell) , 
and an asymptotic flat FLRW region. These models are 
capable of producing good fits to the existing supernova 
data, for an observer at the centre of symmetry. 



IV. RESULTS: SPHERICAL COLLAPSE 

In this section, we evaluate our 4 measures of accel- 
eration in the spherical collapse model (as described in 
Section III A[ ) . This model is known to be able to have an 
accelerating Buchert average, despite locally decelerating 
everywhere [7j. 

In Figure [2] we plot the magnitude of sources, /i, that 
an observer in the spherical collapse model would see, mi- 
nus the magnitude a source at the same redshift would 
have in a pure vacuum model (with all of the dust regions 
removed). The model used in this figure is one in which 
the vacuum regions are 80 Mpc wide at the present time, 
and dust regions are 15 Mpc wide. We refer to these plots 
as 'the Hubble diagram'. We consider two different ob- 
servers in this plot: one at the centre of a collapsing dust 
region (red line) , and one at the centre of an expanding 
vacuum region (blue line). It can be seen that the type 
of region that the observer finds themselves in can have 




Figure 2: Magnitude of sources in the spherical collapse 
model, minus the magnitude they would have in pure vacuum. 
Observations made by a single observer, from the centre of the 
dust and vacuum regions, are displayed as solid red (lower at 
2=1) and blue lines, respectively. The dashed black line is 
the same quantity in the Buchert averaged model. 



a considerable impact on the behaviour of the Hubble 
diagram they construct at low z. At high z, however, 
the curves evolve almost identically (up to a vertical dis- 
placement). In both cases, the blueshifting that occurs 
in the collapsing region forces the curves upwards, due 
to the magnitude continuing to increase even while the 
redshift decreases. This leads to an overall positive gra- 
dient for the curves in Figure [21 while the gradient at any 
point on each individual curve is negative. This result is 
significant, as accelerating FLRW models have positive 
gradient in these plots, while decelerating FLRW models 
have negative gradient. 

In Figure |3] we plot the results of averaging the Hubble 
diagrams constructed by observers in both the dust and 
vacuum regions of the same model considered in Figure [2] 
We do this by constructing individual distance-redshift 
relations for a large number of observers (all in the same 
model), binning these relations in redshift, and then cal- 
culating the mean and standard deviation in each bin. 
We also calculate the magnitudes that would be found in 
the Buchert averaged model, and the spatially averaged 
deceleration parameter ge — 9KS- In the upper panel 
of Figure [3] we plot the magnitude of sources, minus the 
magnitude they would have at the same redshift if all the 
dust regions were removed (as in Figure pi). In the lower 
panel we plot the magnitude of sources minus the mag- 
nitude they would have in the Buchert averaged model. 
In this figure we shift the curves so that they coincide at 
z = 0. This corresponds to a change in the local Hubble 
rate. 

It can be seen from the upper panel of Figure [3] that 
the Buchert average (dashed black line) closely traces 
the mean observed magnitude (solid red line). The 1 
and 2cr confidence regions appear to oscillate, because of 
the jaggedness of the individual distance-redshift curves 
that were averaged over (see Figure [2]). Nevertheless, 
the curve for the Buchert average stays within the la 
confidence region for the entire redshift range considered. 
Part of the reason for this is that the black line has been 
shifted vertically to match the red line at high z, in order 
to aid comparison. Without this shift, the two curves 
have the same zero-point, but fluctuations in the mean 
observational curve at low z cause an offset to build up at 
higher z. The zero-points are the same because they are 
governed by the same spatial average of the local Hubble 
rate {Hq) (for the observational curve) and the Buchert 
average Hubble rate H-p (for the Buchert average curve). 
These are equal at z = 0. 

It can be seen from the lower panel in Figure [3] that 
making the comoving size of the regions smaller, so that 
the initial region is less dominant, significantly reduces 
the fluctuations of the mean observational curve (in both 
cases, these fluctuations decrease as redshift increases). 
This in turn reduces the offset that develops between the 
observational curve and the Buchert average curve. 

The short blue line in the upper panel of Figure [3] is 
the 'effective' curve obtained by setting qq = {qq) in the 
FLRW series expansion for dAiz). Again, qq = qks in the 
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Figure 3: Upper Panel: Magnitude of sources in the spiierical collapse model, minus the magnitude they would have in pure 
vacuum. The solid red line is the mean observed magnitude, obtained by spatially averaging distance-redshift relations on a 
surface of t =constant. The pink bands are la and 2a confidence regions. The black dashed line shows the corresponding 
quantity in the Buchert averaged model. The short blue line is the distance modulus for a FLRW model with deceleration 
parameter q = {qe}- Lower Panel: A difference plot of the red line and dashed line, from the panel above. The solid black line 
corresponds to a model with regions that are a quarter of the size of those used in the upper panel. The pink (wide) and grey 
(narrow) bands are la confidence regions for the two different models. 



spherical collapse model, so this line is also the effective 
curve for the Kristian-Sachs deceleration parameter. It 
can be seen to bear little resemblance to the observational 
and Buchert average curves. 

The values for the deceleration parameters measured 
in this model are 



(qe) 0^ 0.017 
qv ^ -0.167, 



(54) 
(55) 



with (qKs) = {le)j and (gobs) — 1v- By considering other 
model parameters, we have confirmed that the curves due 
to averaging observations, and the Buchert average dis- 
tance curve, always seem to show similar functional be- 
haviour (and thus have similar deceleration parameters) . 
Neither of these quantities are ever close to (gxs) or {qe), 
unless the dust regions are made to expand, or are made 
to be very small. 

We find that the existence of acceleration in both the 
Buchert average and the observational distance-redshift 
relations is a generic feature of spherical collapse models 
consisting of expanding and collapsing regions. Figure |4] 
shows the results of varying the parameters of the FLRW 
regions that constitute the model. The base model has 
the parameters hj = —2.0 and i^mj = 1.8 in the dust 
regions, and hu = 0.7 and flm.n = in the vacuum 
regions. The comoving sizes of these regions are taken 
to be xi = 15 Mpc and xii = 80 Mpc, respectively. 



The collapsing region is chosen to be 52% of its max- 
imum age. The plots in Figure [4] are for models with 
these parameter values, unless they are the parameter 
being varied. The figure shows a strong dependence of 
the Buchert average acceleration on region size - mod- 
els with collapsing regions that are relatively larger have 
greater acceleration, as expected. Increasing the Hubble 
rate in the vacuum regions reduces the current age of the 
model, and correspondingly increases (qs)- This is due 
to this quantity being evaluated at an earlier time, when 
the collapsing regions are more dominant. Figure |4] shows 
that the deceleration parameter of the Buchert averaged 
model is negative as long as the collapsing region takes 
up a non- negligible fraction of the comoving volume. If 
the vacuum region dominates, both the local volume and 
Buchert average deceleration parameters tend to zero. 
Since the observational deceleration parameter is well- 
approximated by the Buchert average measure, it seems 
that all that is required for observers to see an apparent 
acceleration is the existence of a non-negligible fraction 
of collapsing regions. 



V. RESULTS: KASNER-EDS 

We now repeat the analysis of the previous section for 
the Kasner-EdS model. This is included in order to al- 
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Figure 4: Dependence of the Buchert average deceleration parameter, q-v (daslied black line), and spatially-averaged local volume 
deceleration parameter, (qe) (solid red line), on the parameters of the spherical collapse model. Shown is the dependence on 
the Hubble rates in the collapsing and vacuum regions, hi and hu, and the comoving region sizes, xi a-nd xn- 



lay concerns that the coherence in the accelerations of 
the Buchert average and observational measures seen in 
the spherical collapse model were caused simply by inad- 
equacy of the model. 

In Figure [5] we plot the magnitude of sources, /x, that 
an observer in the Kasner-EdS model would see, minus 
the magnitude a source at the same redshift would have 
in a pure vacuum model (with all of the dust regions 
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Figure 5: The same quantities that were plotted in Figure [2] 
but using the Kasner-EdS model. 



removed). The model used in this figure again has dust 
regions that are 15 Mpc wide, but the vacuum regions are 
now ^ 1500 Mpc across. This large number is required in 
order to get cosmologically interesting redshifts, because 
of the rapid acceleration in the vacuum regions along the 
preferred line of sight^. We also choose the model such 
that the dust regions collapse to a singularity in ^ 5 bil- 
lion years. We consider two different observers in Figure 
p] one at the centre of a collapsing dust region (green 
line) , and one at the centre of an expanding vacuum re- 
gion (blue line). Again, it is apparent that the type of 
region that the observer finds themselves in can have a 
considerable impact on the behaviour of the Hubble dia- 
gram they construct at low z, while at high z the curves 
evolve almost identically (up to a vertical displacement). 
Blueshifting again occurs in the collapsing regions, forc- 
ing the curves upward. 

In Figure [6] we plot the results of averaging the 
distance-redshift relations of observers in both the dust 



^ It should be reiterated here that we are not proposing this model 
as a realistic representation of the actual Universe, but rather 
to provide a simple exact solution of Einstein's equations that 
allows us to study relativistic behaviour that may occur in more 
realistic solutions. In this sense, it should be considered as proof- 
of-concept only. 
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Figure 6: The same quantities that were plotted in Figure [3l but using the Kasner-EdS model. 



and vacuum regions of the same model considered in Fig- 
ure [51 as was done in Figure |3] for the spherical collapse 
model. Magnitudes in the Buchert averaged model, and 
the spatially averaged deceleration parameter qq = gKS, 
are also calculated. With this model, however, no verti- 
cal shift of the Buchert averaged model curve is required 
in order for it to agree well with the average of the ob- 
served magnitudes. The Buchert average closely traces 
the mean observed magnitudes. We do, however, perform 
a vertical shift on all of the curves in Figure [6] simultane- 
ously, so that they approach the origin at 2 = 0. Again, 
reducing the comoving size of the regions, so that the 
initial region is less dominant, significantly reduces the 
fluctuations of the mean observational curve. 

The short blue line in the upper panel of Figure [6] is 
again the 'effective' distance modulus curve obtained by 
setting qo = (ge) in the FLRW series expansion for d/i(z). 
We also again have qe — gxs- The values for the decel- 
eration parameters measured in this model are 



qv 



-3.96 

-7.18, 



(56) 

(57) 



with (^Ks) = ile), and (gobs) - Qv- This value of (qe) 
corresponds to rapid acceleration, but not as rapid as 
simply looking through pure vacuum regions in the di- 
rection of inhomogeneity, where qq — —4. The value of 
qx>, on the other hand, corresponds to considerably more 
acceleration than simply looking through the vacuum re- 
gions alone. As in the spherical collapse models, there- 
fore, the presence of the collapsing dust regions causes a 
dramatic increase in both the acceleration of the Buchert 
averaged model, and the average of the observational ac- 



celeration. 

By considering models with other parameter values, 
we again confirm that the curves due to averaging obser- 
vations, and the Buchert average distance curve, always 
seem to show similar functional behaviour (and thus have 
similar deceleration parameters). Neither of these quan- 
tities are ever close to (^ks) or {qe), unless the dust re- 
gions are removed, or become vanishingly small. 



VI. RESULTS: LTB 

The final model that we consider is an LTB space- 
time, with parameters chosen to produce a good fit to 
the supernova data for an observer at r = 0. These are 
displayed in Figure [Tj Unlike the other models we have 
considered, the LTB model has smooth density and Hub- 
ble rate profiles on spatial slices. There are, however, no 
collapsing regions at t = Iq. 

Figure [9] displays the magnitudes of sources in the LTB 
model for our various measures of acceleration. Since 
the chosen model has no homogeneity scale, we plot the 
distance modulus curves for two different averaging do- 
main radii: rp = 1000 Mpc (well inside the void), and 
rx) = 3000 Mpc (near the void boundary) . The results in 
the two cases differ significantly, with the smaller aver- 
aging domain displaying acceleration for the mean obser- 
vational curve (solid red line), while its counterpart for 
the larger domain (solid blue line) shows a strong decel- 
eration. Note that the observational curves are for the 
monopoles of the distance-redshift relation, obtained us- 
ing the dipole approximation described in Section |III C| 
The curves for the Buchert average (dashed lines) do not 
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Figure 7: Density and expansion rates as a function of radius in an LTB void model. The model has spatial curvature of the 
form given in Eq. (53 1, with A^ = —3.82 x 10~*Mpc~^ and Wk = 1800 Mpc. Shown are the density, p (solid black), transverse 



Hubble rate. Hi (blue dashed), and radial Hubble rate, H2 (blue dash-dot). All of the curves are normalised to their values in 
the asymptotically-homogeneous region (p = 9.2 x lO^°M0Mpc~^, H\ — H2 = 57.7 kms~^Mpc~^) 
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Figure 8: Distance modulus curves for an off-centre observer in the same model, at r = 915 Mpc. The dotted blue lines show 
the distance modulus as a function of redshift for radial lines of sight looking out of and into the void (upper and lower curves 
respectively). The dashed blue line is the distance modulus for the monopole of the distance-redshift relation, calculated using 
the dipole approximation. The solid black line is the relation for an observer at the centre of symmetry. 



match the mean observational curves (solid lines) in any 
discernible way. In particular, both of the curves for the 
Buchert average are decelerating, while the averaged ob- 
served relation for the smaller domain size is accelerating. 

The dotted and dash-dotted lines are 'effective' dis- 
tance moduli for the local volume and Kristian-Sachs de- 



celeration parameters. These correspond to series expan- 
sions of the FLRW dA{z) relation with deceleration pa- 
rameters (70 — {<!&) and go — (^Ks) respectively. Neither 
follows the Buchert average, nor the mean observational 
curves. They are closely related to one another though, 
with only a slight discrepancy apparent at higher red- 
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Figure 9: Distance moduli in the LTB model. The red curves (uppermost) are for averaging domains of size r-p — 1000 Mpc, 
and the blue curves for r-p = 3000 Mpc. The solid lines are the mean distance moduli found by spatially averaging the monopole 
of the observed distance-redshift relation, with accompanying la and 2a confidence bands. The dashed lines are the distance 
moduli for the distance-redshift relation in the Buchert averaged spacetime. The dotted and dot-dashed lines are the distance 
moduli for a series expansion of the FLRW dA{z) relation, with deceleration parameters q = (qe) and q = (qkb), respectively. 
The curves have all had their zero-points shifted to match at 2 = 0. 
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Figure 10: Spatial averages of various deceleration parameters as a function of averaging domain radius, r-p. The local volume 
(dash-dotted green line), Kristian-Sachs (dotted blue), and Buchert average (dashed black) deceleration parameters have rather 
similar values. A close-up of these curves is shown in the inset so that they can be distinguished. The thick, solid yellow line 
is the deceleration parameter found by fitting an FLRW model to the distance-redshift relation for a wide range of redshifts, 
z <\. The narrower solid red line shows the same measure, but for a fit to only z < 0.1. 
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shift for the larger domain size. This is because the two 
measures only differ by a term proportional to the shear, 
which is relatively small in this case. 

The upper set of observational curves in Figure |9] show 
that a substantial fraction of observers within a radius of 
r < 1000 Mpc will infer acceleration from the monopole 
of their Hubble diagram, despite all of the other deceler- 
ation parameters indicating deceleration. This is caused 
by spatial variations in expansion rate along their lines of 
sight, but does not seem to be linked with spatial average 
properties of the spacetime, as it was in the spherical col- 
lapse and Kasner-EdS models. Fig. [8]suggests that many 
of these observers will see a large dipole in the distance- 
redshift relation over their sky. Presumably, most of the 
observers would conclude that they lived in a Universe 
that is inhomogeneous on large scales, and would there- 
fore not attempt to fit an FLRW model to their obser- 
vations at all. In that case, the question of whether the 
monopole of the observational relation is accelerating or 
not becomes less of an issue. Attempts to summarise 
the (observational) acceleration of the model with one 
number, the monopole gobs, would fail. Attempting to 
construct a homogeneous Buchert average model would 
probably not be seen as a sensible procedure either^. 

The Buchert average curve in Figure [9] is much closer to 
the effective local volume acceleration and Kristian-Sachs 
curves than to the observational curve, although it still 
deviates from both of them substantially. There are a 
number of reasons why this is the case. The first is that 
the model is not statistically homogeneous, and so the 
Buchert average does not represent the 'typical' condi- 
tions that a light ray would experience traveling through 
the spacetime®. Secondly, the backreaction scalar Qp, 
of Eq. (20), is small in this model, because both the 



shear and the variance of the expansion rate are small. 
As such, the Buchert average deceleration parameter qx> , 
defined in Eq. (21 1, is dominated by the same terms 



that appear in the definitions of the Kristian-Sachs and 
local volume measures, most notably the density. Our 
different measures of acceleration do not exactly reduce 
to one another in this case, but they are rather close (as 
can be seen in Figure 10 inset). We also see that in 



the absence of a sensible way of defining a representative 
smooth model (due to the lack of a homogeneity scale) , 
the mean behaviour of light rays should not be expected 
to correspond to the Buchert averaged model. 



In Figure 10 the various (spatially-averaged) decel- 



eration parameters are plotted as a function of averag- 
ing domain radius, rp. The Buchert average, local vol- 
ume, and Kristian-Sachs deceleration parameters track 



one another rather closely over the whole range of domain 
sizes, for the reasons discussed above. The observational 
curves (thick and thin solid lines) have a rather differ- 
ent behaviour. Both are obtained by calculating qq in an 
FLRW model that has been fit to the monopole of the 
distance-redshift relation as a function of observer posi- 
tion^. The resulting position-dependent value of qn{r) is 
then spatially averaged. The thick yellow curve corre- 
sponds to fitting FLRW models to the distance-redshift 
monopole out to z = 1, and the thin red one to z = 0.1. 
The two are very different. In particular, the fit out to 
higher redshifts shows acceleration for small averaging 
domain radii. Figure [H] shows why this is the case: At 
low redshift, the observational distance modulus curve 
is decelerating (which is also suggested by the fact that 
Qks > everywhere), but appears to accelerate at higher 
redshifts for observers inside the void. By fitting FLRW 
models to the larger redshift range, more of the apparent 
acceleration is captured. This behaviour is wholly due 
to quantities being integrated along past null directions, 
and is not caused by the local curvature of spacetime at 
any one point. 

The results we have presented for the observational 
deceleration are subject to the validity of the dipole ap- 
proximation, that was discussed in Section |IIIC[ This 
approximation will tend to overestimate the monopole of 
the angular diameter distance at a given redshift, causing 
a decrease in gobs relative to its exact value. As can be 
seen from Figure [TOJ this only serves to worsen the dis- 
crepancy between (gobs) and the other deceleration pa- 
rameters at high redshift. At low redshift the approxima- 
tion is most accurate [76 , and so the correction required 
in this regime should be expected to be small. As such, 
we consider our results to be robust to the use of this 
approximation. 



VII. DISCUSSION 

In this paper, we have studied different concepts of 
what it means for a spacetime to display 'accelerating 
expansion'. The measures of acceleration that are asso- 
ciated with these concepts all reduce to the same quan- 
tity in a perfectly homogeneous and isotropic FLRW uni- 
verse, but in an inhomogeneous universe we have shown 
that they can be very different indeed. This occurs even 
to the extent that some can indicate deceleration, while 
others indicate that exactly the same spacetime is accel- 
erating. In universes that are statistically homogeneous 
on large scales, we find, that in order to estimate the 
acceleration inferred by making observations over large 



"■ Instead, the Buchert averaging procedure could be used to de- 
fine a 'smoothed-out' model that is inhomogeneous. This would 
result in position-dependent average quantities, and would likely 
be sensitive to both the shape and size of the averaging domain. 

® In fact, it cannot, as the null geodesies involved in calculating 
this quantity go outside of the averaging domain. 



^ The best-fit FLRW model may not always be a particularly good 
fit - the monopole of the distance-redshift relation in LTB models 
can take on much more complicated functional forms than are 
allowed in FLRW. 
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distances (as is the most usual way to infer acceleration 
in cosmology), one is best off using a model constructed 
from non-local averages of geometric quantities, as oc- 
curs in Buchert's formalism, rather than considering the 
local expansion rate of space. This is in agreement with 
an argument put forward by Rasanen [89l |90j . The ap- 
pearance of acceleration in observations made over large 
scales does not necessarily imply or require the expan- 
sion of space to be accelerating, nor does it require local 
observables to indicate acceleration. 

The models that we used to reach these conclusions are 
given in Section [lll[ and include both exact and approx- 
imate relativistic toy models that are known to display 
some of the types of acceleration we have considered. To 
be specific, we consider: (A) an approximate 'spherical 
collapse' model (with disjoint collapsing and expanding 
FLRW regions); (B) an exact Kasner-EdS model (with 
expanding and collapsing regions along the line of sight); 
and (C) an exact LTB model (expanding everywhere). 
The spherical collapse model has the advantage of be- 
ing able to model reasonably complicated distributions of 
matter, while the Kasner-EdS model allows one to model 
a universe that is statistically homogeneous along the line 
of sight. 

In the spherical collapse and Kasner-EdS models, 
the reconstructed distance-redshift relation, which cor- 
responds most closely to what is actually measured in 
observational cosmology, is closely related to the Buchert 
average, and not the mean local properties of the space- 
time. This means that showing that local spacetime can- 
not accelerate without A, or a quintessence field, is not 
sufficient to disprove backreaction as a source for the ap- 
parent late-time accelerating expansion of the Universe 
[71IH]- This does not, of course, mean that the observed 
acceleration can currently be said to be due to backreac- 
tion: the situations we considered here are very much toy 
models (albeit ones that we expect to capture some of the 
properties of the real Universe). More study is required 
and, in particular, more realistic, non-perturbative mod- 
els of the Universe are required, before any definite con- 
clusions can be drawn about the real Universe. One 
recent attempt at constructing such models appears to 
show some evidence of the effects we describe here [5T] , 
but it is still neither conclusive nor fully realistic. 

As a corollary of our study, a possible source of ob- 
servational evidence for the hypothesis that the appar- 
ent acceleration of the Universe is due to inhomogeneity 
presents itself: If the parameters of FLRW models in- 
ferred from local observations are significantly different 
from those inferred from observations made over large 
distances, then this would seem to imply that the FLRW 



model that we use to model local spacetime is different to 
the FLRW model that best describes the evolution of the 
Universe on large scales. Any such difference would signal 
a significant departure from the predictions of the concor- 
dance ACDM model of the Universe, and would therefore 
cast considerable doubt on the detection of A ^ 0. Of 
course, inferring cosmological parameters from observa- 
tions made on small scales is a considerable challenge. 
Sample variance due to the presence of local structures, 
and the peculiar velocities they induce, would have to be 
very well-understood. Nevertheless, the supposed detec- 
tion of the 'Hubble Bubble' [52] suggests that it may not 
be entirely impossible. 

The link we have found between observations and the 
spatial average can be explained by considering that, for 
a large enough collection of null rays, the typical condi- 
tions experienced by a ray at a given time, t, are likely to 
correspond to the average of local conditions on a hyper- 
surface of t ^constant. These averages are exactly what 
Buchert's approach is constructed to estimate. As long 
as spacetime is statistically homogeneous and isotropic 
above some scale, the result then follows (assuming ob- 
servers and sources are distributed in a volume- weighted 
way). These issues have been considered in detail in |38j . 



In the case of the LTB model, in Section VI we found 
that the Buchert average was more closely related to the 
local measures than to the typical distance-redshift rela- 
tion of an observer, apparently in contradiction with our 
previous results. This model, however, is not statistically 
homogeneous or isotropic, and shows apparent (observa- 
tional) acceleration only for a limited range of averaging 
domain sizes. Furthermore, the distance-redshift rela- 
tion is only poorly represented by its monopole alone, 
since the dipole of the relation would certainly be im- 
portant too. Therefore, our conclusions are that in this 
case the Buchert average is not enough to characterise 
the 'typical' properties of the spacetime, and the mean 
of the monopole of the distance-redshift relation is also 
not enough to characterise what a typical observer should 
expect to see. This leads one to question what it really 
means for an LTB model to exhibit 'average acceleration' 
at all. 
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